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PREFACE 

The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing Is an 8-year program of research, development, evaluation, and 
application of aerospace remote sensing for agricultural resources, which 
began In fiscal year 1980. This program Is a cooperative effort of the 
National Aeronautics and Space Administration, the U.S. Agency for 
Internationa! Development, and the U.S. Departments of Agriculture, Commerce, 
and the Interior. 

The work which Is the subject of this document was performed within the Earth 
Resources Research Division, Space and Life Sciences Directorate, at the 
Lyndon B. Johnson Space Center, National Aeronautics and Space Administration. 
Under Contract NAS 9-15800, personnel of Lockheed Engineering and Management 
Services Company, Inc., performed the tasks which contributed to the 
completion of this research. 
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1 . SCOPE 


This report documents the software which has been developed under Job Order 
73-315, Crop Stage Development Estimation Research. The software consists of 
programs which implement the basic Robertson crop phenology model and the 
Doraiswan^/Thompson crop phenology model. 


2. SYSTEM OVERVIEW 


The main program routines for the Do ralswaiqy /Thompson cr(^ phenology model and 
the basic Robertson crop phenology model are DTMAIN and BRMAIN. The flow 
diagrams for these models are shown In figures 1 and 2. Routines DTMAIN and 
BRMAIN read meteorological data files and coefficient files, accept the 
planting date Information and other Information from the user, and Initiate 
processing. Subroutines DTMOD and BRMOD serve as the driver routines for the 
processing and output segments of the programs. Dally processing for the basic 
Robertson program consists only of calculation of the basic Robertson Increment 
of crop development. Additional processing In the Dora1swan\y/Thompson program 
Includes the calculation of a moisture stress Index and correction of the basic 
Increment of development. Output for both programs consists of listings of the 
dally results. 

All of the routines In these two programs, with the exception of DTMAIN and 
BRMAIN, should be compatible with FORTRAN conpllers on other computer systems. 
The main programs contain FORTRAN- IV plus "OPEN" and "CLOSE" statements. 

The DoralswaiiQr/Thompson crop phenology program and the Robertson crop phenol- 
ogy program are described fully In sections 3 and 4, respectively. Refer to 
Appendix A for descriptions of terms used In this document. Appendix B 
consists of source listings for the entire Dora Iswamy /Thompson program and 
Appendix C consists of source listings for two subroutines unique to the basic 
Robertson program. 
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Figure 2-1.- Flow liagram of the Dora iswamy /Thompson program. 

The numbers showi below the subroutines refer to the section 
numbers in which the subroutines are described in this document. 
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Figure 2-2.- Flow diagram of the Robertson Crop Phenology progrmn. 
The numbers shown below the subroutines refer to the section 
numbers in which the subroutines are described In this document. 
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3. DTMAIM 


The OTMAIN program Is the driver routine of the Dora Iswan^/Thonp son crop 
phenology program. 

CALLING SEQUENCE: PROGRAM DTMAIN 

CALLING ARGUMENTS : Not applicable. 

FILES : Not applicable. 

COMMON BLOCKS : The following abbreviations are used in tables throughout this 
document. 

I - Integer 
R - real 


The common block for the DTMAIN routine is LUN(LUNCOE,LUNRES,LUNMET,LUNTRM, 
LUNRPT) . 


Name 

LUNCOE 

LUNRES 

LUNMET 

LUNTRM 

LUNRPT 


Type Definition 

I Logical unit number of coefficient file. 

I Logical unit number of result files. 

I Logical unit niniber of weather data file. 

I Logical unit number of user terminal. 

I Logical unit number of result report file. 


SUBROUTINES CALLED: 


Name 

Definition 

Reference 

(section) 

METRDR 

Input weather data 

3.1 

INCOEF 

Input coefficients 

3.2 

DTMOD 

Ooraiswamy /Thompson model 

3.3 
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CALLED BY: Not applicable 


LOCAL VARIABLES; 

Name Definition 

COEF Coefficients; C0EF(I,J,K), wnere I equals crop transition period, J 
oquals function (V1,V2,V3), and K equals position of coefficient in 
function; C0EF(*,3,1) is not used. 

CRIT Critical levels; CRIT(I,J), where I equals crop transition period and 
J equals function. 

WD Weather data; WD(I,J), where I equals maximum temperature (**F), minimum 
temperature r*F), and precipitation (cm), and J equals Julian day of 
year. 

OUT Result array; 0UT(I,J), where I equals Julian day of year and J equals 
data field, for example temperature or precipitation. 

NAME File name. 

SM Maximum second layer moisture content. 

PLOTFL Result save flag; If flag equals Y, save; If flag equals N, do not 
save. 

PRTFL6 Print flag; if flag equals Y, print dally results; If flag equals N, 
do not print. 

PLT Julian planting date. 

SEG Segment number. 

LAT Segment latitude. 

AL Constant. 

KOUNT NunRier of days In result array. 
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3.1 HETRDR 


The METRDR subroutine Inputs meteorological data for the current segment. 
CALLING SEQUENCE ; SUBROUTINE METROR(WO^LAT,SEG) 

CALLING ARGUMENTS: 


Name Usage 

Type 

Definition 

WD Output 

R 

Weather data; (1) maximum temperature. **F, (2) minimum 
temperature, *F, (3) precipitation. 

LAT Output 

R 

Segment latitude. 

SEG Output 

R 

Segment ninnber. 

FILES: 

File 


Usage Record fom»t 

Weather data 

file 

Input Variable 


COWPN BLOCKS : Refer to program DTMAIN (section 3) for the description of 

common block LUN. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : The METRDR routine is called by subroutine DTMAIN (section 3). 
LOCAL VARIABLES: Local variables for subroutine METRDR are as follows. 


Name 

Definition 

M 

DO loop index. 

N 

DO 1 oop i ndex . 

YR 

Year. 

IDA 

Flag. 

J 

DO loop index. 
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Name Definition 


JD Julian date. 

ND Julian date. 

K DO loop Index. 
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3 .2 INCOEF 


The INCOEF subroutine Inputs coefficients and calculates critical levels. 

CALLING SEQUENCE: SUBROUTINE INCOEF(COEF.CRIT) 

CALLING ARGUMENTS : 

Name Usage Type Definition 
COEF Output R Coefficients. 

CRIT Output R Critical levels. 

FILES : 

File Usage Record format 
Coefficient Input 8E10.4 

COMMON BLOCKS : Refer to the DTMAIN routine (section 3) for the description of 

common block LUN. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY ; The INCOEF subroutine Is called by subroutine DTMAIN (section 3). 

LOCAL VARIABLES : Local variables for the INCOEF subroutine are as follows. 
Name Definition 

JS DO loop Index; crop stage. 

K DO loop Index. 

L DO loop Index. 
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3.3 OTMOD 


The DTMOD subroutine Is a driver routine for the Dorai swaniy /Thompson crop 
phenology model. Beginning on January 1 and ending when the crop reaches its 
final stage of (tevelopment, the following actions are performed dally. 

1. Water budget and increment of development calculations 

2. Output of daily results if flag is set 

3. Storage of daily results if flag is set 

CALLING SEQUENCE : SUBROUTINE DTMOO(COEF,CRIT,WD,PLT,LAT,AL,SM,PLOTFL,PRTFL6, 

0UT,K0UNT) 

CALLING ARGUMENTS: 


Name 

Usage 

Type 

Definition 

COEF 

Input 

R 

Coefficients . 

CRIT 

Input 

R 

Critical levels. 

WD 

Input 

R 

Weather data. 

PLT 

Input 

I 

Julian planting date. 

LAT 

Input 

R 

Lati tude . 

AL 

Input 

R 

Constant required in transpiration calculations. 

SM 

Input 

R 

Maximum second-layer moisture content. 

PLOTFL 

Input 

R 

Result save flag; if flag equals Y, save; if flag 
equals N, do not save. 

PRTFLG 

Input 

R 

Result print flag; if equals Y, print; if equals N, do 
not print. 

OUT 

Output 

R 

Result storage array. 

KOUNT 

Output 

I 

Number of days in result storage array. 

FILES: 

Not applicable. 
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C0W10N BLOCKS : Refer to main routine DTMAIN (section 3) for the description 
of common block LUN. 

Results for potential (nonlimiting moisture) conditions and observed (limiting 
moisture) conditions are kept separate. Variables holding these results begin 
or end in P (potential) or 0 (observed) and are contained In the unnamed 
common block. 


Name 

Type 

Definition 

JD 

I 

Current Julian date. 

TMIN 

R 

Minimum temperature, ®F. 

TMAX 

R 

Maximum temperature, ®F 

PRECIP 

R 

Precipitation. 

FLB 

R 

(P and 0); first-layer budget. 

SLB 

R 

(P and 0); second-layer budget. 

DAYS 

R 

(.’ and 0); days after stage 1 evaporation. 

K 

R 

(F and 0); transition period. 

CUM 

R 

(P and 0); cumulative increment of development 

LAI 

R 

(P and 0); leaf area index. 

SUMEO 

R 

(P and 0); cumulative bare soil evaporation. 

SUMTR 

R 

(P and 0); cumulative transpiration. 

SUMP 

R 

(P and 0); cumulative precipitation. 

R 

R 

(P and 0); stress ratio. 

SUBROUTINES CALLED: 

Name 


Reference 

Definition (section) 

DTDAY 

Calculations for current day 3.3.1 

PRTDAY 

Print 

results for current day 3.3.2 
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C^LED BY ; The DTMOD subroutine Is called by the DTMAIM driver routine 
(section 3). 

LOCAL VARIABLES ; 

Name Definition 

JDO Julian day when crop reaches final transition period. 

SWITCH Flag used In computing JDO. 
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3 .3 .1 DTDAY 


The DTDAY subroutine calculates the water budget and the Increments of 
development for the current day. In addition, cumulative Increments of 
development and other running totals are updated. 

CALLING SEQUENCE : SUBROUTINE DTDAY( COEF,CRIT,PLT,LAT,AL,SM) 

CALLING ARGUMENTS: 


Name 

Usage 

Type 

Definition 

COEF 

Input 

R 

Coefficients. 

CRIT 

Input 

R 

Critical levels. 

PLT 

Input 

I 

Julian planting date. 

LAT 

Input 

R 

Lati tude . 

AL 

Input 

R 

Constant required In transpiration calculations 

SM 

Input 

R 

Maximum second-layer moisture content. 


FILES : Not applicable. 

COMMON BLOCKS : Refer to the D1M)D subroutine (section 3.3) for the 

description of the common blocks. 

SUBROUTINES CALLED: 


Reference 


Name 

Definition 

(section) 

CORECT 

Adjust dally Increment of 
development for water stress 

3 .3 .1 .5 

OAYLG 

Computes day length 

3 .3 .1 .3 

LAICAL 

Computes leaf area Index 

3 .3 .1 .1 

ROB 

Computes basic Robertson Increment 
of development 

3.3 .1 .4 

WATBUD 

Performs water budget calculations 

3.3 .1 .2 
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CALLED BY : The DTDAY subroutine is called by the DTMOD subroutine 
(section 3.3). 


LOCAL VARIABLES ; 

Name Definition 

DUMl Dumn^ argument . 

DUM2 Dummy argument . 

DUM3 Dummy argument . 

RD Difference between water-limiting and nonlimiting water budget ratios; 

stress index. 

DIDO Uncorrected dally Increment of development; old stage. 

DIDN Uncorrected dally Increment of development; new stage. 

DIDCl Corrected dally increment; old stage. 

DIDC2 Corrected dally Increment; new stage. 

DIDI Dally Increment Interpolated from variables DIDCl and DIDC2. 

DIFF The portion of DIDCl contributing to DIDI when stage crossover occurs. 

RPCT Percentage of DIDC2 contributing to DIDI when stage crossover occurs. 

TMPCUM Temporary cumulative Increment of development used to determine If 
DIDCl carries crop to next transition period. 

DL Day length. 

I TEMP Temporary transition period variable. 

IFLAG Flag used for signaling when subroutine WATBUD should update running 
totals; if equals 0, update; if equals 1. do not update. 

KS Soil -moisture-limiting factor. 
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3.3.1 .1 LAICAL 


The LAICAL subroutine calculates the leaf area Index. 
CALLING SEQUENCE: SUBROUTINE LAICAL(CUM.AL,LAI) 

CALLING ARGUMENTS: 


Name 

Usage 

Type 


Definition 

CUM 

Input 

R 

Cumul ati ve 

increment of development. 

AL 

Input 

R 

Maximum leaf area index; constant for a particular crop 

LAI 

Output 

R 

Cal cul ated 

leaf area index. 


FILES ; Not applicable. 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : The LAICAL subroutine is called by subroutine DTDAY 
(section 3.3.1) . 

LOCAL VARIABLES : 

Name Definition 

BL Constant which is a function of AL (see calling arguments). 
DEV Amount of development since last stage crossover. 

CL Constant which is a function of AL (see calling arguments). 
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3 .3 .1.2 WATBUD 


The WATBUD subroutine performs water budget calculations. On each entry, the 
following actions are performed: 

a. calculation of bare soil evaporation 

b. update of first-layer budget 

c. update of second-layer budget 

d. update of cumulative precipitation 

e. update of cumulative bare soil evaporation 

The following actions are performed only after the crop has been planted: 

a. calculation of transpiration 

b. calculation of stress ratio 

c. update of cumulative transpiration 

CALLING SEQUENCE : SUBROUTINE WATBUD(JD,PLT,TMIN,TMAX,PRECIP,IFLAG,LAI ,AL,SM, 

KS,FLB,SLB,TJ,SUMEO,SUMTR, SUMP, RATIO) 

CALLING ARGUMENTS: 


Name 

Usage 

Type 

Definition 

JD 

Input 

I 

Current Julian date. 

PLT 

Input 

I 

Julian planting date. 

TMIN 

Input 

R 

Minimum temperature, ®F. 

TMAX 

Input 

R 

Maximum temperature, *F, 

PRECIP 

Input 

R 

Precipitation. 

I FLAG 

Input 

I 

If = 1, do not update cumulative totals; 
if = 0, update cumulative totals. 

LAI 

Input 

R 

Leaf area Index. 

AL 

Input 

R 

Constant required In transpiration calculation. 

SM 

Input 

R 

Maximum second-layer moisture content. 

KS 

Input/output 

R 

Soil -moisture-limiting factor (0.0 - 1.0). 
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Name 


Usage Type 


Definition 


FLB 

Input/output 

SLB 

Input/output 

TJ 

Input/output 

SUMEO 

Input/output 

SUMTR 

Input/output 

SUMP 

Input/output 

RATIO 

Output 

FILES: 



File Usage 

Dally results Output 


R First-layer budget. 

R Second-layer budget. 

R Number of days since stage 1 evaporation. 

R Cumulative bare soil evaporation. 

R Cumulative transpiration. 

R Cumulative precipitation. 

R Stress ratio. 


Record format 

Intermediate results for debugging purposes. 


COMMON BLOCKS ; Refer to driver routine DTMAIN (section 3) for a description 
of the common block LUN. 


SUBROUTINES CALLED ; 

Reference 

Name Definition (section) 

PETF Computes potential evaporation 3.3.1 .2.1 

CALLED BY : The WATBUD subroutine is called by the OTDAY subroutine 
(section 3.3.1). 

LOCAL VARIABLES : 

Name Definition 

BLK Leaf area index corresponding to 50 percent of the crop cover, A 
function of AL (see calling arguments). 

ALFA Constant for a particular crop and climatic condition, 

C Constant dependent on soil hydraulic properties. 
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Name Definition 

FLBT First-layer budget for today. This is an estimate of today's first- 
layer moisture content. It is used to determine which second-layer 
budget equation applies. Variable FLBT » FLBY - EO + PRECIP. 

FLBY First-layer budget for yesterday. This is an estimate of yejterdqy's 
first-layer moisture content. It is used in the calculation of the 
FLBT. FLBY is not allowed to drop below zero. 

TAU Temporary variable in bare soil evaporation and transpiration 
calculations . 

TR Transpiration. 

ALFAV Constant for a particular crop. 

ETB Available second-layer moisture; Independent variable In the 
calculation of KS. 

PET Potential evaporation. 

TI Temporary variable In stage-2 bare soil evaporation calculation. 

EOP Temporary variable In stage-2 bare soil evaporation calculation. 

AU Maximum first-layer moisture content. 

EO Bare soil evaporation. 

BL Function of AL (see calling arguments); used In transpiration 

calculation. 


3.3.1 .2.1 PETF 


The PETF function computes potential evaporation. 

CALLING SEQUENCE : FUNCTION PETF(TMIN,TMAX) 

CALLING ARGUMENTS ; 

Name Usage Type Definition 

TMIN Input R Minimum temperature, *F. 

TMAX Input R Maximum temperature, ®F. 

FILES : Not applicable. 

COMMON BLOCKS ; Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY ; The PETF function is called by the WATBUD subroutine 
(section 3.3.1 .2) . 

LOCAL VARIABLES : 

Name Definition 

ES Saturation vapor pressure derived from daily minimum temperature. 
E Actual vapor pressure derived from daily maximum temperature. 
PANEV Pan evaporation. 
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3 .3 .1.3 DAYLG 


The DAYLG function computes day length. 

CALLING SEQUENCE; FUNCTION OAYLG(W.J) 

CALLING ARGUMENTS ; 

Name Usage Type Definition 
W I i{>ui R Latitude. 

J Input I Julian date. 

FILES : Not applicable. 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : The DAYLG function Is called by subroutine OTDAY (section 3.3. 
LOCAL VARIABLES: Not applicable. 
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3.3.1 .4 ROB 


The ROB subroutine performs the dally Increment of (tevelopment calculations 
for the basic Robertson model . 

CALLING SEQUENCE ; SUBROUTINE ROb(TMIN,TMAX,DL.COEF,CRIT,JS,DIDO.OIDN,OIDI.CUH) 
CALLING ARGUMENTS: 


Name 

Usage 

Type 

Definition 

TMIN 

Input 

R 

Minimum temperature, *F. 

TMAX 

Input 

R 

Maximum temperature, “F. 

DL 

Input 

R 

Day length. 

COEF 

Input 

R 

Coefficients. 

CRIT 

Input 

R 

Critical levels. 

JS 

Input/output 

1 

Transition period. Incremented In this subroutine 
If development for today moves crop to next 
transition period. 

DIDO 

Output 

R 

Daily Increment of development using coefficients 
and critical levels from old transition period. 

DIDN 

Output 

R 

Dally increment of development using coefficients 
and criticul levels from new transition period. 

DIDI 

Output 

R 

Dally Increment of development that is Interpolated 
from variables DIDO and DIDN. 

CUM 

Input/output 

R 

Cumulative increment of development. Variable DIDO 
or DIDI (if stage crossover occurs) Is added to CUM 
in '.?vis subroutine. 

FILES 

: Not applicable. 


COMMON BLOCKS: Not applicable. 
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SUBROUTINES CALLED: 


Reference 

Name Definition (section) 

ROBF Computes basic Robertson dally 3.3.1 .4.1 
Increment of development. 

CALLED BY : The ROBF subroutine Is called by subroutine DTDAY (section 3.3.1). 
LOCAL VARIABLES ; 

Name Definition 

TMPCUM Temporar 7 cumulative Increment of development. Used to check if DIDO 
carries crop to new transition period. 

DIFF The portion of DIDO contributing to DIDI when stage crossover occurs. 
RPCT Percentage of DIDN contributing to DIDI when stage crossover occurs. 
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3.3.1 .4.1 ROBF 

The ROBF function computes the basic Robertson dally Increment of development. 
CALL IMG SEQUENCE: FUNCTION ROBF(DL,TMIN,TMAX,COEF,CRIT.JS) 

CALLING ARGUMENTS: 


Name 

Usage 

Type 

Definition 

DL 

Input 

R 

Day length. 

TMIN 

Input 

R 

Minimum temperature, “F. 

TMAX 

Input 

R 

Maximum temperature, ®F, 

COEF 

Input 

R 

Coefficients. 

CRIT 

Input 

R 

Critical levels. 

JS 

Input 

I 

Current transition period. 


FILES: Not applicable. 

COMMON BLOCKS: Not applicable. 

SUBROUTINES CALLED: 

Reference 

Name Definition ( section) 

QUAD Evaluates quadratic. 3. 3. 1.4. 1.1 

CALLED BY: The ROBF function is called by the LAICAL subroutine 

(section 3. 3. 1.1). 
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LOCAL VARIABLES: 


Name Definl tion 

VI Value of quadratic in day length. 

V2 Value of quadratic In maximum temperature. 

V3 Value of quadratic In minimum temperature. 
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3 .3 .1.4 .1.1 QUAD 


The QUAD function evaluates a quadratic. 


CALLING SEQUENCE; FUNCTION QUAD(X,AD,A1,A2,C) 


CALLING ARGUMENTS: 


Name Usage Type 

X Input R 

AO Input R 

A1 Input R 

A2 Input R 

C Input R 


Definition 

Day length or temperature. 
Coefficient. 

Coefficient . 

Leading coefficient. 
Critical level . 


FILES : Not appl i cabl e . 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : The QUAD function is called by subroutine ROBF (section 3.3. 1.4.1). 
LOCAL VARIABLES: Not applicable. 
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3 .3 .1.5 CORECT 


The CORECT subroutine makes an adjustment to the basic Robertson dally 
Increment of development. The corrected dally Increment of development Is a 
function of the current transition period and the moisture stress Index, as 
follows. 

Transition 

period DIDC 

1 DID 

2 DID 

3 DID*EXP(-R) 

4 DID*EXP(R) 

5 DID*F(R) If R < 0.8, or 
010*1.9 if R > 0.8 

where 

R = stress ratio 

DID = daily increment of development before correction 
DIDC = corrected daily increment of development 
F = quadratic function (see source listing) 

CALLING SEQUENCE : SUBROUTINE CORECT(DID,RD,CUM,DIDC) 

CALLING ARGUMENTS ; 

Name Usage Type Definition 

DID Input R Uncorrected daily increment of development. 

RD Input R Stress ratio. 

CUM Input R Cumulative increment of development. 

DIDC Output R Corrected daily increment of development. 

FILES : Not applicable. 
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COMMON BLOCKS ; Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY ; The CORECT subroutine is called by the DTDAY subroutine 
(section 3.3.1) . 

LOCAL VARIABLES ; 

Name Definition 
A1 Constant. 

A2 Constant. 

A4 Temporary variable. 
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3.3.2 PRTDAY 

The PRTDAY subroutine outputs the current day's results for both limiting and 
nonllmiting moisture conditions. 


CALLING SEQUENCE ; SUBROUTINE PRTDAY 
CALLING ARGUMENTS; Not applicable. 


FILES ; 

File Usage Record format 

Dally results Output Two records per day are printed. Record 1 consists of 

Julian date and meteorological data followed by results 
for nonlimiting moisture conditions. Record 2 consists 
only of results for limiting moisture conditions. 

COMMON BLOCKS : Refer to the DTMOD subroutine {section 3.3) for the description 

of common blocks, and refer to the main routine DTMAIN (section 3) for a 
description of common block LUN. 


SUBROUTINES CALLED: Not applicable. 


CALLED BY: The PRTDAY subroutine Is called by subroutine DTMOD (section 3.3). 


LOCAL VARIABLES; Not applicable. 


4. BRMAIN 


The BRMAIN routine Is the main program of the basic Robertson crop phenology 
model. The BRMAIN program Inputs meteorological data (>€TRDR), coefficients 
(INCOEF), and crop planting dates. 

CALLING SEQUENCE; PROGRAM BRMAIN 

CALLING ARGUMENTS; Not applicable. 

FILES; 

File Usage Record format 

Meteorological data Input Refer to METRDR source listing 
Coefficient Input 8E10.4. 

COMMON BLOCKS : Refer to the DTMAIN (section 3) for the description of common 

block LUN. 

SUBROUTINES CALLED: 


Name 

Definition 


Reference 

(section) 

METRDR 

Input meteorological 

data 

4.1 

INCOEF 

Input coefficients 


4.2 

BRMOO 

Basic Robertson model 

driver 

4.3 


CALLED BY: Not applicable. 


4-1 


LOCAL VARIABLES: 


Variable 

Definition 

COEF 

Coefficients; C0EF(I,J,K), where I equals crop transition perioa, J 
equals function (V1,V2,V3), and K equals position of coefficient in 
function; C0EF(*,3,I) is not used. 

CRIT 

Critical levels; CRIT{I,J), where I equals crop transition period 
and J equals function. 

WD 

Weather data; UD(I,J), where I equals maximum temperature (**F), 
minimum temperature (^F), and precipitation (cm), and J equals 
Julian day of year. 

NAME 

File name. 

PLT 

Julian planting date. 

SE6 

Segment number. 

LAT 

Segmenc latitJde. 
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4.1 METRDR 

The METRDR subroutine Is described In section 3.1 of this doctment. 


4-3 


4.2 IKCOEF 

The INCOEF subroutine Is described In section 3.2 of this document. 
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4 .3 BRMOD 


The BRMOD subroutine Is the driver routine for the basic Robertson crop 
phenology model . 

CALIIHG SEQUENCE : SUBROUTINE BRMOO(COEF,CRIT,WD,PLT.LAT) 

CALLING ARGUMENTS: 


Name 

Usage 

Type 

Definition 

COEF 

Input 

R 

Coefficients . 

CRIT 

Input 

R 

Critical levels. 

WO 

Input 

R 

Weather data. 

PLT 

Input 

I 

Julian planting date 

LAT 

Input 

R 

Latitude. 


FILES : 

File Usage Record format 

Report of results Output Each record consists of Julian date, minimum and 

maximum temperature, day length, transition period, 
dally increments of development (old stage, new 
stage, and interpolated), and cumulative increment 
of development. 

COMMON BLOCKS : Refer to the DTMAIN subroutine (section 3) for a description 

of common block LUN. 

SUBROUTINES CALLED : 

Reference 

Name Definition (section) 

DAYLG Computes day length. 3.3.1 .3 

ROB Computes the basic Robertson 3.3.1 .4 
daily increment of development. 


4-5 


CALLED BY : The BWDD subroutine Is called by driver routine BRMAIN 
(section 4) . 

LOCAL VARIABLES : 

Name Definition 


JS 

Transition period. 


TMAX 

Maximum temperature, "F. 


TMIN 

Minimum temperature, *F. 


DL 

Day length. 


DIDO 

Dally Increment of development; 

old stage. 

DIDN 

Dally Increment of development; 

new stik.je. 

DIDI 

Dally increment of development; 
stage. 

interpolated from old stage and new 

CUM 

Cymulative increment of development. 
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4.3.1 DAYLG 


The OAYLG subroutine Is iescribed In section 3.3.1 .3 of this document. 
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4.3.2 ROB 


The ROB subroutine Is described in section 3.3.1 .4 of this document. 
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4.3.2.1 ROBF 

The ROBF subroutine Is described In section 3.3.1 .4.1 of this document. 


4 .3 .2 .1.1 QUAD 

The QUAD subroutine Is described In section 3.3 .1.4 .1.1 of this document. 
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APPENDIX A 


DESCRIPTION OF TERMS USED IN DOCUMENT 


a. Daily Increment of Development (DID) — a numeric representation of a 
crop's development. The DID ranges from 1. at planting to 6. at ripening. 

b. Cumulative Increment of Development (CUM) — the sum of the DID's from 
planting day to current day. 


c. Crop Stage — a point in the crop's development which is identified by one 
of the following events; 


Planting (P) 
Emergence (E) 
Jointing (J) 
Heading (H) 
Soft dough (S) 
Ripening (R) 


(CUM = 1.0) 
(CUM = 2.0) 
(CUM * 3.0) 
(CUM = 4.0) 
(CUM = 5.0) 
(CUM = 6.0) 


e. Transition period _ one of the five intervals between each crop stage, as 
follows. 


p 

to 

E 

(1. 

< 

CUM 

< 

2.) 

E 

to 

J 

(2. 

< 

CUM 

< 

3.) 

J 

to 

H 

(3. 

< 

CUM 

< 

4.) 

H 

to 

s 

(4. 

< 

CUM 

< 

5.) 

S 

to 

R 

(5. 

< 

CUM 

< 

6.) 


f. Stage crossover _ the advancement of a crop to a new transition period. 
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rORTAAN IV-PLUS V02-S1E 
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C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0001 


MAIN PROGRAM 

THIS PROGRAM IS A DRIVER FOR THE OORAISWAMY/THOMPSON 
CROP PHENOLOGY MODEL. 


COMMON BLOCK LUN: 

NAME MEANING 

LUNCOE LOGICAL UNIT NO. OF COEFFICIENT FILE 

LUNRES LOGICAL UNIT NO. OF RESULT FILES 

LUNMET LOGICAL UNIT NO. OF WEATHER DATA FILE 

LUNTRM LOGICAL UNIT NO. OF USER TERMINAL 

LUNRPT LOGICAL UNIT NO. OF RESULT REPORT FILE 


SUBROUTINES CALLED: METRDR , I NCOEF , OTMOD 


LOCAL VARIABLES: 

NAME MEANING 

COEF COEFFICIENTS; COEF(I.J.K) 

I - CROP TRANSITION PERIOD 
J - FUNCTION 

K - POSITION OF COEFFICIENT IN FUNCTION; 
COEF(*.3,l ) IS NOT USED. 

CRIT CRITICAL LEVELS; CRIT<I,0> 

I - CROP TRANSITION PERIOD 
J - FUNCTION 

WO WEATHER DATA; WD(I,J) 

I ‘ FIELD 

1. MAXIMUM TEMPERATURE, DEG. F 

2. MINIMUM TEMPERATURE, DEG. F 

3. PRECIPITATION 
J - JULIAN DAY OF YEAR 

OUT RESULT ARRAY; OUT(I,J) 

I - JULIAN DAY OF YEAR 
J - FIELD 
NAME FILE NAME 

SM MAXIMUM SECOND LAYER MOISTURE CONTENT 

PLOTFL RESULT SAVE FLAG; IF*’ Y', SAVE; IF«’N',DO NOT SAVE. 

PRTFLG PRINT FLAG; IF=‘Y’, PRINT DAILY RESULTS; 

IF»‘N‘ , DO NOT PRINT. 

PLT JULIAN PLANTING DATE 

SEG SEGMENT NUMBER 

LAT SEGMENT LATITUDE 

AL CONSTANT 

KOUNT NUMBER OF DAYS IN RESULT ARRAY 


COMMON /LUN/ L UNCOE , L UNRE S , L UNMET , L UNTRM , L UNRPT 
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FORTRAN IV-PLUS V02-51E 
DTMAIN.FTN /TR:BLOCKS/WR 


0002 

C 

DIMENSION CCEK(5,3,3 >,CRIT(5,3),WD(3,3i 


r 

1 OUT<200,I4> 

0003 

V 

INTEGER PLT.SEG 

0004 

r 

REAL LAT 


C 

INITIALIZATION 

0005 

L. 

LUNCOE - 1 

0006 


LUNRES » 2 

0007 


LUNMET - 3 

0008 


LUNRPT « 4 

0009 


LUNTRM « 5 

0010 


AL « 3.5 

0011 

r 

NAME(30) « 0 


C 

c 

INPUT WEATHER DATA 

0012 

L 

10 

WRITE (LUNTRM, 20) 

0013 

20 

FORMAT (‘STYPE NAME OF WEATHER DATA FILE 

0014 


READ (LUNTRM, 40) ( NAME ( I ) , I * 1 , 29 ) 

0015 

40 

FORMAT (29A1 ) 

0016 


OPEN ( UN I T«L UNMET, NAME -NAME, TYPE *= ’OLD' ,EI 

0017 


CALL METRDR (WD.LAT.SEG) 

0018 


CLOSE (UNIT=LUNMET) 


L 

C 

r- 

INPUT PLANTING DAY 

0019 

L. 

50 

WRITE (LUNTRM, 60) 

0020 

60 

FORMAT (‘STYPE PLANTING DAY (OULIAN) >') 

0021 


READ (LUNTRM, 00) PLT 

0022 

80 

FORMAT ( 14 ) 

0023 


IF (PLT.LT.l .OR. PLT. GT. 366) GO TO 50 


0024 

0025 

0026 

0027 

0028 
0029 


0030 

0031 

0032 

0033 

0034 


0035 


C 

c 

c 

c 

100 


c 

c 

c 

200 

210 

220 

C 

c 

c 

225 


INPUT COEFFICIENTS, CALCULATE CRITICAL LEVELS 
WRITE (LUNTPvM,100) 

FORMAT (‘STYPE NAME OF DORAI SWAMY/THOMPSON COEFFICIENT FILE >') 
PvEAD (LUNTRM,40) < NAME< I ) , I * 1 , 29 ) 

OPEN (UNIT«=LUNCOE,NAME«NAME,TYPE»‘OLD' ,ERR-1000) 

CALL INCOEF (COEF,CRIT) 

CLOSE ( UNIT=LUNCOE ) 

INPUT RESULT FILE INSTRUCTIONS 

WRITE (LUNTRM,2I0) 

FORMAT {‘SGENERATE RESULT FILE ? (Y OR N> >') 

READ (LUNTRM.220) PLOTFL 
FORMAT ( A1 ) 

IF (PLOTFL. NE. 'Y' .AND. PLOTFL . NE N ' > GO TO 200 
INPUT PRINTING INSTRUCTIONS 
WRITE (LUNTRM,230> 
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FORTRAN IV-PLUS V02-S1E 
DTMAIN.FTN /TR :BLOCKS/WR 


8T036 

0037 

0038 


0039 

0040 

0041 

0042 

0043 

0044 


004 S 

0046 

0047 

0048 


0049 


0050 

0051 


230 FORMAT (‘SPRINT DAILY RESULTS 7 (V OR N) >‘) 

READ (LUNTRM,220) PRTFLG 

IF (PRTFLG.NE. • V .AND. PRTFLG . NE . ' N ‘ ) GO TO 225 
C 
C 

C INPUT NAME OF PRINT FILE 

C 

IF (PRTFLG. EQ. 'N* ) GO TO 450 
WRITE (LUNTRM.400) 

400 FORMAT ('STYPE NAME OF PRINT FILE >*> 

READ (LUNTRM,40) ( NAME ( 1 ) , I « 1 , 29 ) 

OPEN (UNIT-LUNRPT.NAME-NAME .TYPE- ' NEW , RECORDS I ZE - 1 33 , ERR- 1000 > 
450 CONTINUE 
C 

C INPUT SECOND LAYER CAPACITY 

C 

WRITE (LUNTRM.460) 

460 FORMAT ('STYPE SECOND LAYER CAPACITY >') 

READ (LUNTRM.480) SM 
480 FORMAT ( F 10. 1 > 

C 

C 

c 

c 

CALL DTMOD ( COE F , CR I T . WD , PLT . L AT , AL , SM , P L OTF L . PRTF LG , 

1 OUT.KOUNT) 

C 

c 

c 

c 

IF (PRTFLG. EO. 'V ) CLOSE ( UN I T*L UNRPT ) 

C 

IF (PLOTFL.EQ. 'N* ) GO TO 499 
C 

C CREATE FILE OF RESULTS 

C 


0052 


DO 490 0*2. U 

0053 


OPEN (UNIT-LUNRES,NAME*'DTPL0T.DAT' ,TYPE-'NEW ) 

0054 


WRITE (LUNRES.485 ) ( OUT( I . 1 ) . I » 1 , KOUNT ) 

0055 


WRITE (LUNRES.485 ) ( OUT( I , 0 ) . 1 » 1 . KOUNT ) 

0056 

485 

FORMAT (F10.Z) 

0057 


CL?SE ( UNIT-LUNRES > 

0058 

490 

CONi ^NUE 

0059 

L 

KOUNT * KOUM * 2 

0060 


WRITE ( LUNTRM.491 ) KOUNT 

0061 

491 

FORMAT (•0'.I3.* RECORDS IN EACH RESULT FILE.*) 


c 

c 


006? 

C 

499 

WRITE (LUWTRM.500) 

0063 

500 

FORMAT (’SCONTINUE 7 ( Y OR N ) > ' ) 

0064 


READ (LUNTRM.220) FLAG 

006 b 


IF (FLAG.NE.'V .AND. FLAG.NE.'N') GO TO 499 

0066 


IF (FLAG.EQ, ‘ Y’ ) GO TO 10 

0067 


STOP 
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FORTRAN IV-PLUS V02-51E 
DTMAIN.FTN /TR : BLOCKS/WR 

C 

C 

C 


0068 

1000 

CONTINUE 

0069 


WRITE (LUNTRM,102ier) 

0070 

1020 

FORMAT (‘TERROR IN OPENING 

0071 


STOP 

0072 


END 
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FOmAN IV-PLUS V02-51E 
METRDR.FTN /TR : BLOCKS/WR 


0001 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE METRDR( WO , LAT , SEC ) 

THIS SUBROUTINE INPUTS METEOROLOGICAL DATA FOR CURRENT 
SEGMENT. 


CALLING ARGUMENTS: 


NAME 

I/O 

MEANING 

WD 

0 

WEATHER DATA 

1. MAXIMUM TEMPERATURE 

2. MINIMUM TEMPERATURE 

3. PRECIPITATION 

LAT 

0 

SEGMENT LATITUDE 

SEG 

0 

SEGMENT NUMBER 


CALLED BY: DTMAI N , BRMA I N 


LOCAL VARIABLES: 

NAME 

M 

N 

0 

K 

YR 

IDA 

JD 

ND 


MEANING 

DO LOOP INDEX 
DO LOOP INDEX 
DO LOOP INDEX 
DO LOOP INDEX 
YEAR 
FLAG 

OULIAN DATE 
JULIAN DATE 


ROB01580 


0002 


COMMON /LUN/ LUNCOE , L UNRE S . L UNMET , L UNTRM , L UNRPT 


0003 

c 

DIMENSION W0( 3,366 ),ME r( 16 > 

ROB01590 

0004 


INTEGER T( 2 . 19 ) ,SEG. YR 

ROB01600 

0005 


REAL LAT 

ROB016I0 

0006 


DO 1 N « I , 366 

ROB01620 

0007 


DO 1 K « 1 .3 

ROB01630 

000 0 


WD(M.N) * 9999. 

RO&01640 

0009 

1 

CONTINUE 

ROB01650 

0010 


READ( LUNMET.9 ) SEG. LAT 


0011 

9 

FORMAT( I4,11X.F5.2> 

ROB01670 

0012 


I » 1 

ROB01680 


C 

C»* INPUT PRECIPITATION 


001 

3 

C 

c 

5 

INPUT IN HUNDRETHS OF AN INCH, CONVERT TO CENTIMETERS 
READILUNMET.gflf.END-ieO) JD.MET.VR 


001 

4 

90 

FORMAT(6X. 13. IX. 16U.4X. 12 ) 

ROC01700 

001 

5 


IF( JD.EQ. -99 )GO TO 20 

ROB01710 
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FORTIUN IV-PLUS V02-51E 
IMCOEF.FTH /TRsBLOCUCS/WR 

fffl SUBROUTINE INCOEF (COEF«CRIT) 

C 

C 

C THIS SUBROUTINE INPUTS COEFFICIENTS AND CALCULATES 

C CRITICAL LEVELS. 

C 

C 

C CALLING ARGUMENTS: 

C 


c 

NAME 

I/O 

MEANING 

c 





c 

COEF 

0 


COEFFICIENTS 

c 

CRIT 

0 


CRITICAL LEVELS 

c 





c 





c 

CALLED BY: DTMAIN 



c 





c 





c 

LOCAL VARIABLES 

: 



c 





c 

NAME 



MEANING 

c 





c 

OS 

DO 

LOOP 

INDEX; CROP STAGE 

c 

K 

DO 

LOOP 

INDEX 

c 

L 

DO 

LOOP 

INDEX 


C 

C 

C 


0002 



COMMON /LUN/ LUNCOE . LUNRES . LUNMET , LUNTRM • LUNRPT 

0003 



DIMENSION COEF(5,3,3),CRIT(5,3) 

0004 



DO 100 0S=1 .5 

0005 



DO 50 K =1,3 

0006 

50 


CRIT(JS,K) = 99999.9 

0007 

L 


READ (LUNCOE,75) ( ( COEF < 0S , K , L ) * L» 1 , 3 ) , K» 1 , 2 ) , COEF(JS,3,2> 



♦ 

COEF< JS,3,3) 

0008 

75 


FORMAT (8E10.4) 

0009 

C 


IF (COEF(OS,1,3).NE.0.0> CRIT(0S,l) » 



♦ 

-<COEF(OS, 1 ,2)~2.0*COEF( JS, 1 , 3 )*COEF ( JS , 1,1)) 



♦ 

/<2.0*COEF<OS,1,3>) 

0010 



IF (COEF(OS,2,3 >.NE.0.0) CRIT(0S,2) = 



♦ 

-(COEF< JS,2,2 >-2.0*COEF( JS,2,3)*COEF(OS,2, 1 ) ) 



♦ 

/( 2.0*COEF(OS,2,3 ) ) 

0011 



IF (COEF<OS,3,3).NE.0.0) CRIT(0S,3> = 



♦ 

-(COEF(OS.3,2 )-2.0*COEF( JS.3,3 )*COEF(OS,2, 1 > > 



4- 

/<2*0*COEF< JS,3,3)) 

0012 

100 


CONTINUE 

0013 



RETURN 

0014 



END 
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OTMOO • FTN /TR : BLOCKS/WR 

0001 SUBROUTINE DTMOD < COEF »CRIT ,W0 , PLT , LAT ,AL , SM . PLOTFL , PRTFLG * 

1 OUT.KOUNT) 

C 

C 

C THIS SUBROUTINE IS A DRIVER FOR THE DORAISWAMY / THOMPSON 

C CROP PHENOLOGY MODEL. STARTING ON JANUARY 1* AND ENDING WHEN 

C THE CROP REACHES ITS FINAL STAGE OF DEVELOPMENT, THE FOLLOWING 

C ACTIONS ARE PERFORMED DAILY: 

C 

C 1. WATER BUDGET AND INCR. OF DEVELOPMENT CALCULATIONS 

C 2. OUTPUT OF DAILY RESULTS, IF FLAG SET 

C 3. STORAGE OF DAILY RESULTS, IF FLAG SET 

C 
C 

C CALLING ARGUMENTS: 

C 


c 

NAMF 

I/O 

MEANING 

L. 

c 

COEF 

I 

COEFFICIENTS 

c 

CRIT 

I 

CRITICAL LEVELS 

c 

WD 

I 

WEATHER DATA 

c 

PLT 

I 

JULIAN PLANTING DATE 

c 

LAT 

I 

LATITUDE 

c 

AL 

I 

CONSTANT REQUIRED IN TRANSPIRATION 

c 



CALCULATIONS 

c 

SM 

I 

MAXIMUM SECOND LAYER MOISTURE CONTENT 

c 

PLOTFL 

I 

RESULT SAVE FLAG. IF«’Y‘, SAVE; IF=’N 

c 



DO NOT SAVE . 

c 

PRTFLG 

I 

RESULT PRINT FLAG. I F» ' Y ‘ , PRI NT ; IF«'| 

c 



DO NOT PRINT. 

c 

OUT 

0 

RESULT STORAGE ARRAY 

c 

KOUNT 

0 

NUMBER OF DAYS IN RESULT ARRAY 


C 

C 

C 

C DESCRIPTION OF UN-NAMED COMMON BLOCK: 

C 

C NOTE: RESULTS FOR POTENTIAL ( NON -L I M I T I NG MOISTURE) CONDITIONS 

C AND OBSERVED (LIMITING MOISTURE) CONDITIONS ARE KEPT SEPARATE. 

C VARIABLES HOLDING THESE RESULTS BEGIN OR END IN P (POTENTIAL) 

C OR 0 (OBSERVED). 

C 


C 

NAME 



MEANING 

c 

C 

JO 

CURRENT 

JULIAN DATE 

c 

TMIN 

MINIMUM 

1 TEMPERATURE, DEG. F 

C 

TMAX 

MAXIMUM 

1 TEMPERATURE. DEG. F 

c 

PRECIP 

PRECIPITATION 

C 

FLB 

( P 

AND 

0); 

FIRST LAYER BUDGET 

c 

SLB 

( P 

AND 

0); 

SECOND LAYER BUDGET 

c 

DAYS 

( P 

AND 

0) ; 

DAYS AFTER STAGE 1 EVAPORATION 

c 

K 

( P 

AND 

0 ) ; 

TRANSITION PERIOD 

c 

CUM 

( P 

AND 

0); 

CUMULATIVE INCR. OF DEVELOPMENT 

c 

LAI 

( P 

AND 

0); 

LEAF AREA INDEX 

c 

SUME0 

( P 

AND 

0); 

CUMULATIVE BARE SOIL EVAPORATION 

c 

SUMTR 

(P 

AND 

0); 

CUMULATIVE TRANSPIRATION 
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C 

C 

R 

<P AND 0>) STRESS RATIO 


c 

c 

c 

SUBROUTINES CALLED: DTOAV, FRIDAY 


c 

c 

c 

CALLED BY: DTMAIN 


c 

c 

LOCAL VARIABLES 



c 

c 

NAME 

MEANING 


c 

JOO 

JULIAN DAY WHEN CROP REACHES 

FINAL TRANSITION 

c 


PERIOD 


c 

SWITCH 

FLAG USED IN COMPUTING JDO 



0002 

r 

COMMON /LUN/ L UNCOE , LUNRES . L UNMET , L UNTRM . LUNRPT 

0002 


COMMON JD,TMIN,TMAX,PRECIP ,PFLB,OFLB.PSLB,OSLB. 


C 

♦ POAYS ,ODAYS ,KP ,KO,PCUM,OCUM,PLAI ,OLAI * 

♦ PSUME0,OSUME0,PSUMTR,OSUMTR ,PSUMP ,OSUMP ,RP,RO 

0T004 

INTEGER PLT 

0005 


REAL LAT 

0006 

r 

INTEGER SWITCH 

0007 

C 

r 

DIMENSION COEF(5,3,3),CRIT(5,3),WD( 3,366 >,OUT(200,U> 


c 

r 

PERFORM INITILIZATIONS AT BEGINNING OF YEAR 

0008 

L 

SWITCH = 0 

0009 


KOUNT = 0 

0010 


PFLB ^ 20.0 

001 1 


0FL3 = 20.0 

0012 


PSLB * SM 

0013 


OSLB * SM 

0014 


PDAYS = 1 .0 

0015 


ODAYS = 1.0 

0016 


PSUMP = 0.0 

0017 


OSUMP = 0.0 

0018 


KP 1 

0015 


KO = 1 

0020 


PCUM = 1 .0 

0021 


OCUM = 1,0 

0022 


PSUME0 = 0.0 

0023 


OSUME0 = 0.0 

0024 


PSUMTR = 0.0 

0025 

r 

OSUMTR = 0.0 

0026 


DO 500 JD=1 ,366 

002 7 


TMAX * W0< 1 , JD ) 

0028 


TMIN = V/D(2,JD) 

0029 


PRECIP « WOO. JO) 


I 
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C 
C 

0030 CALL OTDAY < COEF ,CRIT, PLT , LAT , AL , SM ) 

C 

C 

C PRINT CURRENT DAY’S RESULTS 

C 

0031 IF (OO.GE.PLT .AND. PRTFLG . EQ . ’ Y * ) CALL PRTOAY 
C 

C STORE RESULTS FOR CURRENT DAY 

C 


0032 


IF (PLOTFL.EO. 'W .OR. OD.LT.PLT) 60 TO 400 

0033 


KOUNT - KOUNT ♦ 

1 

0034 


IF ( KOUNT. LE. 200) GO TO 300 

0035 


WRITE (5,280) 


0036 

28f 

FORMAT (• *. ’SUBSCRIPT OF ARRAY, OUT. EXCEEDS LIMIT') 

0037 


STOP 


0038 

300 

CONTINUE 


0039 


OUT( KOUNT, 1) * 

FLOAT(OD) 

0040 


0UT(K0UNT,2) • 

PLAI 

0041 


0UT<K0UNT,3> « 

OLA I 

0042 


0UT(K0UNT,4) * 

PCUM 

0043 


OUTCKOUNT.G) * 

OCUM 

0044 


0UT(K0UNT,6) = 

PSUMP 

0045 


0UT(K0UNT,7) * 

OFLB 

0046 


0UT(K0UNT,8> * 

OSLB 

0047 


0UT<K0UNT,9) = 

OSUME0 

0048 


OUT( KOUNT, 10) = 

OSUMTR 

0049 


OUT( KOUNT, 11 ) = 

PSUME0 

0050 


OUT( KOUNT, 12 ) » 

PSUMTR 

0051 


OUT ( KOUNT, 13 ) » 

RP 

0052 


IF (OCUM.GE.6.0) RO » 0.0 

0053 


OUT( KOUNT, 14) 

RO 

0054 

400 

CONTINUE 


0055 

c 

IF (SWITCH. EQ.l 

.OR. OCUM. LT. 6.0) GO TO 450 

0056 


SWITCH = 1 


0057 


JDO » OD 


0058 

450 

C 

CONTINUE 


0059 

C 

IF (PCUM.GE.6.0 

.AND. OCUM. GE. 6.0) GO TO 600 

0060 

500 

C 

CONTINUE 


0051 

C 

600 

r' 

CONTINUE 


0062 

L 

V/RITE (LUNTRM.800) JDO.JO 

0063 

800 

FORMAT <///'0' , 

•CUM. INCR. OF OEVEL. REACHES 6.0 ON DAY 


♦ 13, • UNDER STRESS CONDITIONS,' 

♦ /’0','CUM. INCR. OF DEVEL. REACHES 6.0 ON DAY 

♦ 13. ‘ UNDER NON-STRESS CONDITIONS,' 

+ /// ) 

C 

C 

0064 RETURN 
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0065 END 
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0001 SUBROUTINE DTOAY < COEF ,CRIT, PLT, LAT , AL , SM ) 


THIS SUBROUTINE PERFORMS WATER BUDGET AND DAILY INCREMENT OF 
DEVELOPMENT CALCULATIONS FOR THE CURRENT DAY. IN ADDITION, 
CUMULATIVE INCREMENTS OF DEVELOPMENT AND OTHER RUNNING TOTALS 
ARE UPDATED. 

BEFORE PLANTING DATE, WATER BUDGET CALCULATIONS ARE THE SAME 
UNDER BOTH WATER-LIMITING AND NON-LIMITING CONDITIONS. 
THEREFORE, THE WATER BUDGET IS CALCULATED ONLY UNDER 
NON-LIMITING CONDITIONS UNTIL PLANTING DAY IS REACHED, 

AT WHICH TIME BOTH BUDGETS ARE EQUATED. ON DAYS AFTER THE 
PLANTING DAY, BUDGETS FOR BOTH CONDITIONS ARE PROCESSED 
SEPARATELY, 

ON ALL DAYS AFTER THE PLANTING DAY, TWO INCREMENTS OF 
DEVELOPMENT ARE COMPUTED. THE FIRST IS THE BASIC 
ROBERTSON INCREMENT. THE SECOND IS DERIVED FROM THE FIRST 
BY MAKING ADJUSTMENTS FOR THE EFFECTS OF WATEPx STRESS. 


CALLING ARGUMENTS: 



NAME 

USAGE 

MEANING 


COEF 

I 

COEFFICIENTS 


CRIT 

I 

CRITICAL LEVELS 


PLT 

I 

JULIAN PLANTING DATE 


LAT 

I 

LATITUDE 


AL 

I 

CONSTANT REQUIRED IN 
CALCULATION 

TRANSPIRATION 

SM 

I 

MAXIMUM SECOND LAVER 

MOISTURE CONTENT 


SUBROUTINES CALLED: LAICAL .WATBUD , DAYLG . ROB ,CORECT 


LOCAL VARIABLES: 

NAME MEANING 


DUMl 

DUM2 

DUM3 

RD 

DIDO 

DIDN 

DIDCl 

DIDC2 

DIDI 

DIFF 

RPCT 


DUMMY VARIABLE 
DUMMY VARIABLE 
DUMMY VARIABLE 

DIFFERENCE BETWEEN WATER-LIMITING AND 
NON-LIMITING WATER BUDGET RATIOS; STRESS INDEX 


UNCORRECTED 

DAILY 

INCR. 

OF 

DEVEL. 

; OLD 

STAGE 

UNCORRtCTED 

DAILY 

INCR. 

OF 

DEVEL. 

; NEW 

STAGE 

CORRECTED 

DAILY 

INCR. 

OF 

DEVEL. 

; OLD 

STAGE 

CORRECTED 

DAILY 

INCR. 

OF 

DEVEL . 

; NEW 

STAGE 


DAILY INCREMENT INTERPOLATED FROM DIOCl 
AND DI0C2 


PORTI 

ON 

OF 

OIOCl 

CONTRI 

BUTING 

TO 

DIOI 

WHE 

STAGE 

CROSSOVER 

OCCURS 





PERCE 

NT 

OF 

0IDC2 

CONTRI 

BUTING 

TO 

DIDI 

WHE 
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C STAGE CROSSOVER OCCURS 

C TMPCUM TEMPORARY CUMULATIVE INCREMENT OF DEVELOPMENT 

C USED TO CHECK IF DIDCl CARRIES CROP TO NEXT 

C TRANSITION PERIOD 

C DL DAY LENGTH 

C ITEMP TEMPORARY TRANSITION PERIOD VARIABLE 

C IFLAG FLAG TO SIGNAL WHEN SUBROUTINE WATBUD SHOULD 

C UPDATE RUNNING TOTALS; IF-0. UPDATE AND IF«1, 

C DO MOT UPDATE 

C KS SOIL MOISTURE LIMITING FACTOR 

C 
C 
C 

0002 INTEGER PLT 

0003 REAL LAT.KS 

0004 DIMENSION COEF ( 5 , 3 , 3 ) , CR I T( 5 , 3 ) 

0005 COMMON JD,TMIN,TMAX,PRECIP .PFLB,OFLB,PSLB,OSLB, 

♦ PDAYS,ODAYS,KP,KO,PCUM,OCUM,PLAI ,OLAI , 

♦ PSUME0.OSUME0J»SUMTR ,OSUMTR,PSUMP *OSUMP ,RP,RO 
C 

C 

C 

e**«*«* BEFORE AND ON PLANTING DAY - WATER BUDGET FOR BOTH LIMITING 

C AND MON-LIMITING CONDITIONS 

C AFTER PLANTING DAY - WATER BUDGET FOR ONLY NON-LIMITING 

C CONDITIONS 

C 

0006 CALL LAICAL ( PCUM, AL ♦ PLAI ) 

0007 KS - 1.0 

0008 IFLAG » 0 

0009 CALL WATBUD ( 0 D , PLT , TM I N , TMAX , P REC I P , I F L AG , P L A I , AL , SM , KS , P F L B , 

♦ PSLB ,PDAYS .PSUME0.PSUMTR ,PSUMP ,RP > 

C 

C 

0010 IF (OD-PLT) 600,200,^00 
C 

C****** ON PLANTING DAY. EQUATE RUNNING TOTALS 

C 

0011 200 OSLB = PSLB 

0012 OFLB = PFLB 

0013 ODAYS « POAYS 

0014 OSUME0 * PSUME0 

0015 OSUMTR « PSUMTR 

0016 OSUMP « PSUMP 

0017 GO TO 600 

(;>*•••• AFTER PLANTING DAY. COMPUTE DAILY INCREMENT OF DEVELOPMENT 

C 

C01Q 400 CONTINUE 

U019 OL = DAYLGI LAT.OD ) 

C 

C*«* MODEL UNDER NON-LIMITING MOISTURE CONDITIONS 

C (WATER BUDGET PROCESSED ABOVE) 

C STAGE ESTIMATE 

C020 CALL ROB ( TM I N . THAX . Dt , COEF ,CR IT , KP . OUMl . DUMB . DUM3 , PCUM ) 

C 
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0021 


0022 

0023 

0024 
002 B 


0026 
002 7 
0028 


0029 

0030 

0031 


C 

c 

c 


c 


c 

c 

CW#* 

c 


/TR:BLOCKS/WR 

MODEL UNDER LIMITING MOISTURE CONDITIONS 
IF (OCUM.GE.6.0) GO TO 600 
WATER BUDGET 

CALL LAICAL ( OCUM ,AL ,OLAI ) 

Kt - 0.0 
IrLAG - 0 

CALL WATBUD ( JD , P LT , TMI N * TMAX , PREC I P . I FLAG , OL AI , AL , SM , KS , 

► OFLB.OSLB ,ODAYS,OSUME0,OSUMTR,OSUMP,RO) 

STAGE ESTIMATE 
ITEMP • KO 
TMPCUM « OCUM 

CALL ROB (TMIN,TMAX,DL,COEF,CRIT, ITEMP, OIOO*DIDN,DIDI, TMPCUM) 

STRESS CORRECTION, OLD STAGE 
RD =* RP - RO 

CALL CORECT ( D I DO , RD , OCUM , D I DC 1 > 

IF (OCUM+OIDCl .GE. FLOAT(KO+l)) GO TO 500 
UPDATE CUMULATIVE D.I.D. 



C 


CURRENT 0. I .D. DUES NOT 

0032 


OCUM 

* OCUM + DIDCl 

0033 


GO 

TO 600 


L 

c 


NEW STAGE ENTERED 

0034 

500 

KO * 

KO + 1 

0035 


TMPCUM » OCUM ♦ DIDI 

0036 


CALL LAICAL ( TMPCUM .AL .OLA I ) 

0037 


KS « 

0.0 

0038 


IFLAG 

; = 1 

0039 


CALL 

WATBUD (OD,PLT.TMIN,TMAX,PRE( 


c 

* 

OFLB.OSLB.ODAVS.DUMl 


L 

c 


STRESS CORRECTION 

0040 


RD « 

RP - RO 

0041 


CALL 

CORECT ( DIDN.RD, TMPCUM, DIDC2 


K. 

C 


INTERPOLATE 

004 2 


DIFF 

= FLOAT( KO ) - OCUM 



RPCT 

« 1 . - Dirf /DIDCl 

0.J4 4 


OCUM 

= OCUM + <DIFK + RPCT**DIDC2> 

0045 

Q 

RETURN 

::/* 4 f; 


Kr:i) 
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SUBROUTINE LAICAL ( CUH,AL .LAI) 

THIS SUBROUTINE CALCULATES LEAF AREA INDEX. 


CALLING ARGUMENTS: 
NAME I/O 


MEANING 

CUMULATIVE INCREMENT OF DEVELOPMENT 
MAXIMUM LEAF AREA INDEX? CONSTANT 
CALCULATED LEAF AREA INDEX 


CALLED BY: DTDAV 


LOCAL VARIABLES: 


MEANING 

- CONSTANT WHICH IS A FUNCTION OF AL(SEE CALLING 
ARGUMENTS > 

- CONSTANT WHICH IS A FUNCTION OF AL 

- AMOUNT OF DEVELOPMENT SINCE LAST STAGE 
CROSSOVER 


0002 


REAL 

LAI 






0003 


LAI - 

0.0 






0004 


BL • 

AL / 

3 

.0 




0005 


IF (C 

UM.LT. 

7 

.0) 

GO 

TO 

500 

0006 


IF (CUM.LT. 

3 

.0) 

GO 

TO 

100 

0007 


IF (CUM.LT. 

4 

.0) 

GO 

TO 

200 

0008 


(CUM.LT. 

5 

.0) 

GO 

TO 

300 

0009 


IF (CUM.LT. 

6 

.0) 

GO 

TO 

400 

0010 


GO TO 

500 






001 I 

100 

DEV • 

CUM *■ 


2.0 




0012 


LAI • 

0.85* 

BL'DEV 



0013 


GO TO 

500 






0014 

200 

DEV - 

CUM - 


3.0 




0015 


LAI ■ 

0.85* 

B 

L^2 

. 15* 

BL- 

^DEV 

0016 


GO TO 

500 






0017 

300 

DEV = 

CUM - 


4.0 




0018 


LAI - 

AL-0. 

6 

•BL’ 

►DEV 



0* T 19 


GO TO 

500 






0020 

400 

DFV « 

CUM - 


5.0 




0021 


CL « - 

AI.-BL* 

0 

.6 




00 22 


LAI « 

CL-CL 

* 

DEV 




0023 

500 

CONTI 

NUE 






hiJcA 


RETUR 

N 






0025 


END 
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0001 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r* 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE WATBUD < OD , PLT.TMIN , TMAX , PREC IP , I FLAG , LAI , AL . SM . 
♦ KS,FLB,SLB,TJ ,SUME0,SUMTR. SUMP , RATIO) 

THIS SUBROUTINE PERFORMS WATER BUDGET CALCULATIONS. 

ON EACH ENTRY, THE FOLLOWING ACTIONS ARE PERFORMED: 

1. CALCULATION OF BARE SOIL EVAPORATION 

2. UPDATE OF FIRST LAYER BUDGET 

3. UPDATE OF SECOND LAYER BUDGET 

4. UPDATE OF CUMULATIVE PRECIPITATION 

5. UPDATE OF CUMULATIVE BARE SOIL EVAPORATION 

THE FOLLOWING ACTIONS ARE PERFORMED ONLY AFTER THE 
CROP HAS BEEN PLANTED: 

1. CALCULATION OF TRANSPIRATION 

2 . CALCULATION OF STRESS RATIO 

3. UPDATE OF CUMULATIVE TRANSPIRATION 


CALLING ARGUMENTS: 


NAME 

I/O 

MEANING 

OD 

I 

CURRENT JULIAN DATE 

PLT 

I 

JULIAN PLANTING DATE 

TMIN 

I 

MINIMUM TEMPERATURE, DEG. F 

TMAX 

I 

MAXIMUM TEMPERATURE, DEG. F 

PRFCIP 

I 

PRECIPITATION 

IFLAG 

I 

IF*1, DO NOT UPDATE SUB-TOTALS; 
IF=^0, UPDATE SUB-TOTALS 

LAI 

I 

LEAF AREA INDEX 

AL 

I 

CONSTANT !N TRANSPIRATION CALCULATION 

SM 

I 

MAXIMUM SECOND LAYER MOISTURE CONTENT 

KS 

I 

SOIL MOISTURE LIMITING FACTOR 

FLB 

I/O 

FIRST LAYER BUDGET 

SLB 

I/O 

SECOND LAYEPv BUDGET 

TJ 

I/O 

DAYS SINCE STAGE 1 EVAPORATION 

SUME0 

I/O 

CUMULATIVE BARE SOIL EVAPORATION 

SUMTR 

I/O 

CUMULATIVE TRANSPIRATION 

SUMP 

I/O 

CUMULATIVE PRECIPITATION 

RATIO 

0 

STRESS RATIO 


SUBROUTINES CALLED: PETF 


CALLED BY: DTDAY 


LOCAL VARIABLES: 

NAME MEANING 

RLK LEAF AREA INDEX CORRESPONDING TO 50. PERCENT 

CROP COVER. A FUNCTION OF AL (SEE CALLING ARGU.) 
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c 

ALFA 

CONSTANT FOR A PARTICULAR CROP CLIMACTIC 

c 


CONDITION 

c 

C 

CONSTANT DEPENDENT ON SOIL HYDRAULIC PROPERTIES 

c 

FLBT 

FIRST LAYER BUDGET FOR TODAY, THIS IS AN 

c 


ESTIMATE OF TODAY ‘ S FIRST LAYER MOISTURE CONTENT. 

c 


IT IS USED TO DECIDE WHICH SECOND LAYER BUDGET 

c 


EQUATION APPLIES. 

c 

TAU 

TEMPORARY VARIABLE IN BARE SOIL EVAPORATION AND 

c 


TRANSPIRATION CALCULATIONS, 

c 

TR 

TRANSPIRATION 

c 

ALFAV 

CONSTANT FOR A PARTICULAR CROP 

c 

ETB 

AVAILABLE SECOND LAYER MOISTURE; INDEPENDENT 

c 


VARIABLE IN CALCULATION OF KS, 

c 

PET 

POTENTIAL EVAPORATION 

c 

TI 

TEMPORARY VARIABLE IN STAGE TWO BARE SOIL 

c 


EVAPORATION CALCULATION, 

c 

E0P 

E0 PRIME; TEMPORARY VARIABLE IN STAGE TWO 

c 


BARE SOIL EVAPORATION CALCULATION. 

c 

AU 

MAXIMUM FIRST LAYER MOISTURE CONTENT 

c 

E0 

BARE SOIL EVAPORATION 

c 

BL 

FUNCTION OF AL (SEE CALLING ARGUMENTS); 

c 

c 


USED IN TRANSPIRATION CALCULATION. 


C 

c 


0002 

n 

COMMON /LUN/ LUNCOE , L UNRES , L UNMET , LUNTRM , LUNRPT 

0003 


REAL KS*LAI 

0004 


INTEGER PIT 

0005 


DATA ALFA,ALFAV/1 . 15, 1 .26/ 

0006 

r 

DATA C,AU/3. 0,20.0/ 

0007 


BL = AL/3.0 

0008 

c 

TR = 0.0 


c 

r 

COMPUTE POTENTIAL EVAPORATION 

0009 

L 

n 

PET « PETF(TMIN,TMAX) 


L 

c 

COMPUTE TRANSPIRATION IF AFTER PLANTING DAY 

0010 

V.# 

TAU =* EXP(~0.537*LAI ) 

001 1 


IF ( JO.LT.PLT+1 1 GO TO 200 


c 

IF PROCESSING WATER BUDGET UNDER NON-LIMITING WATER CONDITIONS 


c 

THEN KS = 1.0 ON ENTRY. COMPUTE KS ONLY WHEN PROCESSING UNDER 


c 

LIMITING WATER CONDITIONS <KS=0.0 ON ENTRY) 

^012 


IF ( KS.EO. 1 .0) GO TO 60 

0013 


ETB = ( SL8/SM)*100.0 

0014 


KS = 1 .0 

0.015 


IF ( ETB.GT.G0.0) GO TO 60 

r/016 


KS = “0.00549764 -0 . 0004 9547*( ETB > +0 . 001 23765 *< < ETB ) >*’^2 ) 
1“0.000022U*( ( ETB )**3 ) ♦ 0.0000001 l*( ( ETB ) 

0017 


IF (KS.GT. 1 .0) KS = 1 .0 

00 1 8 

60 

CONTINUE 

0019 


RLK = BL*2.0 

0020 


IF (LAI.GT.SLK) GO TO 100 
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C 

CROP COVER LESS THAN OR EQUAL TO 50 PERCENT 



0021 


TR = ALFAV*»< 1 .“TAU)*PET*KS 



0022 


GO TO 150 




C 

CROP COVER GREATER THAN 50 PERCENT 



0023 

10J? 

TR = ( 1 .35-TAU )’*'PET*KS 



0024 

150 

IF (TR.GT.PET) TR^PET 



0025 

200 

CONTINUE 




L 

COMPUTE BARE SOIL EVAPORATION 



0026 


FLBY = FLB 



0027 


IF (FLBY. LT. 0.0) FLBY = 0.0 





FLBT = FLBY + PRECIP 



0029 


IF (FLRT.GT.AU) GO TO 300 




L 

c 

/** 

STAGE 2 EQUATION 



0030 

c 

IF ( PRECIP .NE ,0.0) GO TO 280 



0031 


T1 = TJ - 1.0 



0032 


IF (TI.LT.0.0) TI = 0.0 



0033 


GO TO 290 



0034 

280 

E0P = SUME0 - PRECIP 



0035 


IF ( E0P -LT.0.0) E0P = 0.0 



0036 


TI = (E0P/C)**2 



0037 

290 

CONTINUE 



0038 


E0 ^ C*TO**0.5 - C*TI**0.5 



r.03O 


ir ( f0.lt. 0.0) E0 = 0.0 





TT = TJ + 1.0 



'*■ /'» 4 i 


CO TO 350 




L 

C 

p 

STAGE 1 EQUATION 



004 2 

300 

CONTINUE 



004 3 


E0 = TAU’^PET^C 1 .0/ALFA) 



0044 

C 

TJ = 1.0 



0045 

C 

350 

C 

CONTINUE 




C 

c*** 

r 

UPDATE RUNNING TOTALS 




C 

MOTE: PRECIP WHICH IS IN UNITS OF CM ARE ADDED 

TO 


c 

BUDGETS WHICH ARE IN UNITS OF MM. HOWEVER. 

THIS 


c 

OISCREPAMCY IS ACCOUNTED FOR IN WATER BUDGET 

EQUATIONS 

004 6 

c 

IF ( IFLAG.NE .0) GO TO 900 




V. 

C 

SAVE YESTERDAY'S FIRST LAYER BUDGET 



0047 


FLBY = FLB 




c 

FIRST LAYER BUDGET - CUMULATIVE 



00 4 8 


FLB = FLB - E0 + l-HECIP 




c 

FIRST LAVER BUDGET - ACTUAL 




c 

(FIRST LAYFR GUCGET CAN NOT PHYSICALLY BE LESS 

THAN 0. 

0049 


IF (FLBY. LT. 0.0) FLBY = 0.0 



0050 


FL8T = FLBY - E0 + PRECIP 




c 

SECOND LAYER BUDGET 
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0/TB 1 
0053 


0054 

0055 
00 36 


0057 

0058 


0059 

0060 


0061 


:f U3 
'•■05 4 

"1065 


0066 
O0G 7 
006 8 
i !‘!G3 
00 70 

'507 1 


0072 
007 3 


C 

C 

c 

400 


c 

c 

450 



C 

c 


c 


IF (FLBT.CT.AU) GO TO 400 
IF < FLBT .LT.0,0> GO TO 450 

IF (FLBT.GE.0.0 .AND. FLBT.LE.AU) GO TO 500 

FIRST LAYER BEYOND MAXIMUM WATER HOLDING CAPACITY. 

EXCESS WATER ADDED TO SECOND LAYER BUDGET. 

SLB SLB - TR + (FLB-AU) 

FLD = AU 
GO TO 670 

HOT ENOUGH WATER IN TOP LAYER TO SATISFY BARE SOIL EVAPORATION. 
REQUIRED \/ATER IS TAKEN FROM SECOND LAYER. 

SLB = SLB - TR - E0 
GO TO 60.0 

i;0 EXCHANGE DETWEFN FIRST AND SECOND LAYERS. 

THE ONLY SECOHD LAYER LOSS IS TRANSPIRATION. 

SLB » SLB - TR 

CONTINUE 

WATER BEYOND MAX. HOLDING CAPACITY IS LOST AS DRAINAGE 
IF (FLB.CT.SM) SLB = 5M 

UPDATE EVAPORATION. TRANSPIRATION. AND PRECIPITATION TOTALS 
SUMTR ^ SUMTR + TR 

r;iiiTc;V - suriF0 + t:;* 

SUMi' ^ SUMR + PRECIR 

IF { JD.LT.PLT+1 ) GO TO 1000 


C 

COMPUTE RATIO 
C 

900 CONTINUE 

IF {PET.EO.0.0) PET = 0.0000001 
RATIO = (TR + E0) / PET 
IF ( RATIO. LT. 0,0) PkATIO = 0.0 
IF ( RATIO.GT. 1 .0) RATIO = 1.0 
C 

1000 CONTINUE 

r 

C WRITE STATEMENTS FOR DEBUGGING 

C 

D l/RITE ( LUNRPT. 1005 > P F T , TAU . E TB , KS . BL K , TR . E0 . TI 


D1005 

FORMAT ( ‘ 

■ . • PET 

= 

• .Ffi.2.2X. 

D 


' TAU 

= 

' .F6.P.2X, 

D 

+ 

‘ E7B 


’ .F6.2.2X. 

D 

+ 

■ KS 

= 

■ .F6.2.2X, 

D 


• ULK 


' .Ff) .2 .2X . 

D 


‘ TR 


■ ,rr .2.2X, 

0 

♦ 

■ E0 

= 

• .F6.2.2X, 

D 

+ 

• TI 

= 

■ , F G . 2 ) 


PvETURN 

END 


B-27 


FORTRAN IV-PLUS Vff2-ME 10:35'?4 02-JUN-81 PAGE 

WATBUO.FTN / f R : BLOC KS /WP. 


OF POOH * = ^ 


t, ii ?: *: . ■- 

C' CT- S! V - 


.J Js» Q ”3 Cd =■ t: 

cl UJ O — I li. K 


« ft « « 

^ 05 ^ ^ rj 

— — »3 ^ n 'n 

^ c< ■:5 «a Q «B 5J 

Ot ^ S) ^ ^ SS5 ^ 
Bj tl 3« Q C5 S5 S) 
55 ^ ^ 


f\j rsi 'M -t ^ 

ft ft ft ft ft « ft 

□£ Od ^ Q£ Qt 


(J uJ Cli 

< r: Lj 

-1 3K ^ 5: 

>a u- rs -I =» o 

< Ui ^ O- t/l h- 


t: ■S‘ ??> fsi ^ 

(R *<r in "s ro M \o 

<5 5» S{ C3 O Q tj 

tS) &5 5S> *5 Q 

B> B5 BJ ^ B) 5S 

<5 2R B) £5 <5 5:5 (5 


^ -itf fSJ ^ 

ft ft ft ft <- ft ft 

od Qd Qd — Qd ttd a; 


O _> U U <5 

-j _i «j _i eJ 

n: z z 2 q; 
o o o o > 

U O U U O 


!3 5 :5 ::» 

Od cd Qd Od Od 


> O 
< > o 

U. C3 CO 2 ^ 


l5<S}«tvay3&^'«T 

B565&5C&B}B5^5^ 

<S5B5fc?e!555Q50t". 

I I I i 1 t I I 


ftftftftftftftft 

odododod — adodcd 



LU 


UJ 

< K H 

\X3 m 55 m 

S 


£ 

u CO •-• z CD 

Ni n 

< 


< 

-J _J < Z> _J < ^ 

ft*' 

•«♦ 69 VO ^ rj 

z 

655 

z 

<'ou_i_icnh-»- 

« ft ft ft ft *- 

55'«rvD-<t5fS>vDxX5 

fVJ V fV ^ 

{A 

63 

cn 

(VJ-’^Csir-j-'^roQ 

IT 69 55 55 55 

lA 

69 

(A 

69690065696969 

— 69 55 <S} C/ 

UJ 

6S5 

UJ 

69690e5«0690 

55 55 tR 55 55 

Qd 

69 

Qd 

690055^; 0^0 

55 55 65 69 Q 

O 


o 

00055650690 


n 

1 

o 

1 i 1 1 1 1 1 i 

^ < < 

< 


< 

u'oru.uvouuu 


3: o < < cd 

< o o Q < 2 

2 U a. ^ > Z' 

W W W « -1 


fidodfldodi^adfidod 


r ^ rg rn ^ vD 


O. O Qd 
Od H Z 
id efi Z F- r ^ 
_i- 4 «jcn 3 <r 5 Z 
ftC(£)U.id.^OdC0H 


eo <5 GQ 
< «) cr CQ 

—I u> (vj lo 


B-28 


TOTAL SPACE ALLOCATED * 001700 480 



^ ^ ^ 


^ CM r> in u) 


«aN «o4 «<^ 

(M 

(SI (M CM (M CM CM 


^ QB m Qi) Si Q» 

ta 

flO £3 m OQ CO CO 

o 

c o o o o o 


D£ c^ e^ c^ q£ ac 







< 

K -J 

« « 



>- 


LU 

< 3 

ro vo 



< 


r 

-J r> 

«<• «-• 

< 


o 




VD ^ 

CO 






fs. •-* 

• • 


tn 




• • 

in 


UJ 




^ CO 

CO 


H 




♦ 

•• 


=> 


o 


z 

S) 


a. 

• ♦ 

'v, 

*— * 

vn < 



2T 

l/) 



^ H 


r> 

O 

1- 



U) ^ Z « 


m 

(J 

21 



* «< «: r^ 


:» 


Lli 



«— CTt 3 ro 



z 

z 



w » h- * 

> 

c5 

o 

=> 

LJ 


■f CM LU CO 



M 

C5 

2: 


m 1 fii w 


> 

1- 

a: 

< 


»-3 + 


< 

u 

< 

z 


* * • ^ 

u 

o 

z 




CM (M 

LU o 


4> 

C5 



— • ^ • 

^ -J 

z 

U- 

Z 



^ . (M 

in cc 

o 


t-^ 



« Ck — 

1 ♦• 


LO 

—1 



* (JJ 11 z 

CM cr 

H 

»— • 

_J 



ej la • o oi 

<Si H 

u 

X 

< 



^ ^ -J 3 

> 

z 

H 

o 





o 





< O LU < UJ 

(/> 

U- 





O U ^ O OC 







■¥ 








Q. 







> 







Z 


C 

c 

c 

c 

c 

c 

o u 

C 

C 

c 

U 

H 







Z U. 







< • 







CJ ij 







H _i 






CM CO 'O' uT 

c;i >- 






65 isi la o 

O < 






65 55 O is 

L. O 

t; 





6. te ts la 


o 

z 

LU 


VD 


B-29 





1 

[ 


wO 

Ul 

O 

r» 

< 


in 

UJ 

o' 

c. 

t? 

< 


Ul 

r. 


X 

< 


Ul 

X 

< 


f 




CO 

CO 

c/1 

CO 

UJ 

Ul 

tc 


a 

o 

a 

o 

< 

< 


fL 

>■ 


UJ 

X 

< 


Ui 

£ 

< 


CO 

to 

CO 

t/1 

UJ 

UJ 

a: 

0£ 

o 

o 

o 

o 

< 

< 

Ul 

UJ 

cu 

Q. 

>- 

>- 

H 

H 





UJ 

CO 


-J _i _i -I 

£ 

1 


U U L) U 

< 

z 

to 

-J _i — 1 _l 

z 


UJ 

» » • » 



1- 

z z z z 


1 

ID 

o o o o 



CO 

U L> U O 



»— « 

• » > > 

to 


or 

—•003 

CO 


K 

V » * » 

Ul 


H 

> > :5 

Qt: 

*«» 

ro 

< 

Qc; QC or n: 

o 

o 

< 


Ul 

£ 

< 


CO 

UI 

or 

o 

o 

< 


« 

fsj 



FORTRAN 

ROB.FTN 

M£T1 


0002 

0003 

0x;r-4 


IV-PLUS V02-51E 10:35:39 02-JUN-81 PAGE I 

/TR:BLOCKS/WR 

SUBROUTINE ROB ( TMI N ,TMAX , OL .COEF ,CRIT, JS , DI DO, D1 DN , DID I ,CUM ) 
C 

C THIS SUBROUTINE PERFORMS DAILY INCREMENT OF DEVELOPMENT 


c 

c 

CALCULATIONS FOR BASIC 

ROBERTSON MODEL. 

c 

c 

CALLING ARGUMENTS: 


c 

c 

NAME 

I/O 

MEANING 

c 

c 

TMIN 

I 

MINIMUM TEMPERATURE, DEG. F 

c 

TMAX 

I 

MAXIMUM TEMPERATURE, DEG. F 

c 

DL 

I 

DAY LENGTH 

c 

COEF 

I 

COEFFICIENTS 

c 

CRIT 

I 

CRITICAL LEVELS 

c 

JS 

I/O 

TRANSITION PERIOD. INCREMENTED IF 

c 



TODAY’S DEVELOPMENT MOVES CROP TO 

r 



NEXT TRANSITION PERIOD. 

c 

01 DO 

0 

DAILY INCREMENT OF DEVELOPMENT USING 

c 



COEFKICIENTS AND CRITICAL LEVELS FROM 

c 



OLD TRANSITION PERIOD 

c 

01 DM 

0 

DAILY INCREMENT OF DEVELOPMENT USING 

c 



COEFFICIENTS AND CRITICAL LEVELS FROM 

c 



NEW TRANS! ION PERIOD 

c 

c: 

OIDI 

0 

DAILY INCREMENT OF DEVELOPMENT 
INTERPOLATED FROM DIDO AND DIDN 

c 

CUM 

I/O 

CUMULATIVE INCREMENT OF DEVELOPMENT. 

c 



DIDO OR DIDKIF STAGE CROSSOVER OCCURS) 

c 

c 



IS ADDED TO CUM. 

c 

c 

c 

SUBROUTINES CALLED: RORF 

c 

c 

c 

CALLED BY: DTDAY 


c 

c: 

c 

c 

LOCAL VARIABLES 

: 


MAMC 


MEANING 

c 

c 

TMPCUM 

- TEMPORARY CUMULATIVE INCREMENT OF DEVELOPMENT. 

c 


USED 

TO CHECK IF DIDO CARRIES CROP TO NEW 

c 


TRANSITION PERIOD. 

c 

OIFF 

- TH_ PORTION OF DIUO CONTRIBUTING TO DIDI WHEN 

c 


STAGE 

CROSSOVER OCCURS. 

c 

RPCT 

- PERCENT OK DIDN CONTRIBUTING TO DIDI WHEN 

c 

c 


STAGE 

CROSSOVER OCCURS 

c 

DIMENSION COEF( 

5,3,3 ). 

CRIT(5,3) 


C 

C D.I.O. FOR CURRENT STAGE 

DIDO * ROBF( DL ,TMIN,TMAX,COEF ,CRIT,US ) 
THPCUM = CUM + 01 no 
C 
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0005 

0006 

0007 

0008 


0009 

0010 

0011 

0012 


0013 


00U 

0015 

0016 
0/ri7 
0i^ 1 n 
0019 
0/720 
0021 



C 

c 

200 

C 

C 

C 

C 

C 


IF (TMPCUM.LE.6.0) GO TO 150 
CUM ■ 6.0 
OS - OS ♦ 1 
RETURN 


CONTINUE 

IF <TMPCUM.GE.FLOAT(OS+1 )) GO TO 200 

TODAY’S DEVELOPMENT DOES NOT MOVE CROP TO NEXT STAGE 
CUM ■ CUM + DIDO 
RETURN 


CONTINUE 


TODAY'S DEVELOPMENT CARRIES CROP TO NEXT STAGE. 

COMPUTE DAILY INCREMENT OF DEVELOPMENT USING NEXT STAGE 
COEFFICIENTS, THEN INTERPOLATE. 

OS » OS ^ 1 

DIOM « ROBF(DL,TMIN,TMAX,COEF,CRIT,OS) 

DIFF « FL.OAT(OS) - CUM 
RPCT » 1 .0 - DIFF/DIDO 
DIDI = DIFF + RPCT*0I0N 
CUM = CUM + DIOI 
RETURN 
END 
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FORTRAN IV-PIUS V02-51E 
CORECT . FTN /TR : BLOCK5/WR 


0001 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0002 

0003 

0004 

0005 

C 

0006 

0007 

0000 

0009 


SUBROUTINE CORECT ( D ID , RD ,CUM, D I DC ) 


THIS SUBROUTINE MAKES AN ADJUSTMENT TO THE BASIC ROBERTSON 
DAILY INCREMENT OF DEVELOPMENT. THE CORRECTED DAILY 
INCREMENT OF DEVELOPMENT IS A FUNCTION OF THE TRANSITION 
PERIOD IN WHICH THE CROP IS CURRENTLY IN. AS WELL AS THE 
STRESS RATIO. 

TRANSITION PERIOD DIDC 


1 

2 

3 

4 

5 


DIO 

DID 

DID*EXP(-RD) 

DID*EXP(RO) 

DID*C 


WHERE: 

C • F(R) IF R.LE.0.8 AND 1.9 IF R.GT.0.8 
F - QUADRATIC FUNCTION (SEE BELOW) 


CALLING ARGUMENTS: 


NAME 

USAGE 

MEANING 

DID 

I 

UNCORRECTED DAILY INCREMENT OF DEV. 

RD 

I 

STRESS RATIO 

CUM 

I 

CUMULATIVE INCREMENT OF DEVELOPMENT 

DIDC 

0 

CORRECTED DAILY INCREMENT OF DEV. 


CALLED BV: DTDAY 

LOCAL VARIABLES: 

NAME MEANING 


A1 CONSTANT 

A2 CONSTANT 

A4 TEMPORARY VARIABLE 


DATA Al .A2 /I .0. I .0/ 

DIDC • DID 

IF (RD.LE .0.0) GO TO 1000 
IF (CUM. LT. 2.0) GO TO 1000 

IF (CUM.GE .3.0) GO TO 50 
IF (RD.GT.0.25) Al « Al*0.95 
DIDC • DID * ( Al*EXP( -RD ) ) 

GO TO 1000 
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0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 


0019 

0020 


C 

50 CONTINUE 

IF (CUM. CE. 4.0) GO TO 100 
DIDC - DID • (A2*EXP(RD)> 

GO TO 1000 
C 

100 CONTINUE 

IF (CUM. GE. 6.0) GO TO 1000 

A4 ■ (0.9907143 ♦ 0.1869047*RD ♦ 0. 726 1 905*( RD**2 ) ) •• 1.16 
IF (RO.GT.0.8) A4 - 1.9 
DIOC • DID * A4 
C 

C IF(DSTI .LT.0. 1 )A4«l . 

C IF( OSTI .GE .0.0.ANO.OSTI .LT.0.2)A4-l .0 

C IF( DSTI .GE .0.2.ANO. DSTI . LT . 0 . 3 )A4 • 1 . 05 

C IF( DSTI .GE .0.3. AND .DSTI . LT .0 . 4 )A4- I . 1 

C IF( DSTI .GE .0. 4.ANO.DSTI . IT . 0 . 5 )A4 ■ 1 . I 5 

C IF( DSTI .GE .0.5.ANO. DSTI .LT. 0.6 )A4»1 . 2 

C ir( DSTI .GE .0.6. AMD. DSTI . LT . 0 . 7 )A4 - 1 . 35 

C IF( DSTI .GE .0.7.AND.DSTI .LT.0.8 )A4-1 . 45 

C IF(DSTI.GE.0.8)A4«1.6 

C 

1000 RETURN 
END 


PAU02920 

PAU02930 

PAU02940 

PAU02950 

PAU02960 

PAU02970 

PAU029U0 

PAU02990 

PAU03000 
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/TR:BLOCKS/WR 

FUNCTION PETF (TMIN.TMAX) 

C 

C THIS FUNCTION COMPUTES POTENTIAL EVAPORATION. 


C 

C CALLING ARGUMENTS: 


C 

C 

C 

C 


NAME 

I/O 

MEANING 




TMIN 

I 

MINIMUM 

TEMPERATURE, 

DEG. 

F 

TMAX 

I 

MAXIMUM 

TEMPERATURE . 

DEG. 

F 


CALLED BY; WATBUD 


LOCAL VARIABLES; 


NAME 


MEANING 


ES 

E 

PANEV 


SATURATION VAPOR PRESSURE DERIVED FROM 
DAILY MINIMUM TEMPERATURE 
ACTUAL VAPOR PRESSURE DERIVED FROM 
DAILY MAXIMUM TEMPERATURE 
PAN EVAPORATION 


PETF = 0.0 

IF (TMIN.LT.0.0 .OR. TMAX . L E . 32 .0 ) RETURN 
ACTUAL VAPOR PRESSURE (E) 

E =6.11*EXP((-1.762042621E05 ♦ 5.597607915E03*TMAX- 
12.850972636*TMAX**2)/( 1 .2b4162E05 + 273.0*TMAX)) 
SATURATION VAFOR PRESSURE (ES) 
rS = 6.11*EXP((-1.'^62042621E05 ♦5.597607915E03*TM1N- 
12.850972G36’"TMIN**2 )/( I .254162E05 + 273.0*TMIN)) 

PAN EVAPORATION (PANEV) 

PANEV = (0.2163 + 0.3473*E - 0 . 2644 *ES )*( 24 . 55/3 1 .0 > 
POTENTIAL EVAPORATION (PETF) 

PETF = PANEV * .85 

RETURN 

END 
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ROBF.FTN /TR:BLOCKS/WR 


0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
001 1 


FUNCTION ROBF ( JL ,TMIN ,TMAX ,COEF ,CRIT , JS > 

C 

C THIS FUNCTION COMPUTES THE BASIC ROBERTSON DAILY INCREMENT 

C OF DEVELOPMENT. 

C 

C 

C CALLING ARGUMENT: 

C 


C 

c 

C 

NAME 

I/O 

MEANING 


DL 

I 

DAY LENGTH 


c 

TMIN 

I 

MINIMUM TEMPERATURE 

. DEG. F 

c 

TMAX 

I 

MAXIMUM TEMPERATURE 

, DEG. F 

c 

COEF 

I 

COEFFICIENTS 


c 

CRIT 

I 

CRITICAL LEVELS 


c 

c 

OS 

I 

CURRENT TRANSITION 

PERIOD 

c 

c 

c 

SUBROUTINES CALLED: QUAD 


c 

c 

c 

CALLED BY: ROB 




c 

c 

r* 

LOCAL VARIABLES: 




C 

c 

NAME 


MEANING 


c 

VI 

VALUE 

OF QUADRATIC IN DAY 

LENGTH 

c 

V2 

VALUE 

OF QUADRATIC IN MAX. 

TEMPERATURE 

c 

V3 

VALUE 

OF QUADRATIC IN MIN. 

TEMPERATURE 


C 

c 

c 

DIMENSION COEF(5,3,3),CRIT(5,3) 

C 

IF (COEF(OS,1 , 1 ) .NE .999. > GO TO 20 
VI » 1.0 
GO TO 50 


C 

20 VI * QUAD (DL,COEF(OS,1,1),COEF(JS,1,2).COEF(OS,1,3),CRIT<OS,1)> 

50 V2 -QUAO(TMAX,CO£F< JS,2, 1 ) . COEF ( US , 2 , 2 > , COEF ( OS , 2 , 3 ) ,CRIT( JS , 2 ) > 

V3 =QUAD(TMIN,COEFtOS.2, 1 ) , COEF ( OS , 3 , 2 ) . COEF ( OS , 3 , 3 ) ,CRIT( OS , 3 ) ) 
C 

ROBF = VI * (V2 ♦ V3) 


C 


RETURN 

END 
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FORTRAN IV-PLUS, V02-51E 
BRMAIN.FTN /TR : BLOCKS/WR 


0001 

0002 

0003 

0004 


0005 

0006 


C 

C THIS ROUTINE IS THE MAIN PROGRAM OF THE BASIC 

C ROBERTSON CROP PHENOLOGY MODEL. BRMAIN INPUTS 

C METEOROLOGICAL DATA (METRDR), COEFFICIENTS (INCOEF), 

C AND CROP PLANTING DATE, AND THEN RUNS THE MODEL (BRMOD). 

C 

C 

C SUBROUTINES CALLED: METRDR , I NCOEF . BRMOD 

C 

C 

C COMMON BLOCK LUN: 

C 


c 

p 

NAME 

MEANING 




c 

LUNCOE 

LOGICAL UNIT NO, OF 

COEFFICIENT 

FILE 

c 

LUNRES 

LOGICAL UNIT NO. OF 

RESULT FILES 


c 

LUNMET 

LOGICAL UNIT NO. OF 

WEATHER DATA 

FILE 

c 

LUNTRM 

LOGICAL UNIT NO, OF 

USER TERMINAL 


c 

c 

LUNRPT 

LOGICAL UNIT NO. OF 

RESULT REPORT 

FILE 

c 

c 

LOCAL VARIABLES 

: 




c 

p 

NAME 

MEANING 




c 

COEF 

COEFFICIENTS* COEF( 

I ,J,K) 



c 


I - CROP TRANSITION PERIOD 



c 


0 - FUNCTION 




c 


K - POSITION OF 

COEFFICIENT 

IN 

FUNC* 

c 


C0EF(’*,3,1 > 

IS NOT USED. 



c 

CRIT 

CRITICAL LEVELS! CRIT(l.J) 



c 


I - CROP TRANSITION PERIOD 



c 


J - FUNCTION 




c 

WD 

WEATHER DATA; WD< I , 

0 ) 



c 


I - FIELD 




c 


1. MAXIMUM 

TEMPERATURE , 

DEG. F 

c 


2. MINIMUM 

TEMPERATURE, 

DEG. F 

c 


3. PRECIPITATION 



c 


J - JULIAN DAY OF YEAR 



c 

NAME 

FILE NAME 




c 

PLT 

JULIAN PLANTING DATE 



c 

SEG 

SEGMENT NUMBER 




c 

lAT 

SEGMENT LATITUDE 





C 

C 

C 

COMMON /LUN/ L UNCOE , L UNRE S , L UNME T , L UNTRM , L UNRPT 
C 

DIMENSION COEF(5,3,3 ) ,CRIT( 5.3 ) , WD ( 3 , 366 ) . NAME ( 30 ) 
INTEGER PLT.SEG 
REAL LAT 
C 

C INITIALIZE 

C 

LUNCOE « ) 

LUNMET = 3 
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BRMA 1 N . FTN /TR : 8LOCKS/WR 

0007 LUNRPT ■ ^ 

0008 LUNTRM » 5 
C 

C INPUT WEATHtR DATA 

C 

0009 WRITE (LUNTRM, 20) 

0010 20 FORMAT ('STYPE NAME OF WEATHER DATA FILE >') 

0011 READ (LUNTRM, 40) NAME 

0012 40 FORMAT (30A1) 

0013 OPEN (UNIT-LUNMET, NAME-NAME, TYPE •’OLD’ , ERR- 1000) 

0014 CALL METRDR (W0,LAT,SEG) 

0015 CLOSE (UNIT-LUNMET) 

C 

C INPUT PLANTING DAY 

C 

0016 WRITE (LUNTRM, 60) 

0017 60 FORMAT (’STYPE PLANTING DAY (OULIAN) >’) 

0018 READ (LUNTRM, 80) PLT 

0019 80 FORMAT (14) 

C 

C BASIC ROBERTSON MODEL 

C 

0020 WRITE (LUNTRM, 100) 

0021 100 FORMAT ( * STYPE NAME OF BASIC ROBERTSON COEFFICIENT FILE >•) 

0022 READ (LUNTRM, 40) NAME 

0023 OPEN ( UN IT-LUNCOE . NAME-NAME , TYPE • ’ OLD ’ .ERR-1000 ) 

0024 CALL INCOEF (COEF,CRIT) 

0025 CLOSE (UNIT-LUNCOEF ) 

C 

0026 WRITE (LUNTRM, 120) 

0027 120 FORMAT (‘STYPE NAME OF DAILY OUTPUT FILE >‘) 

0028 READ (LUNTRM, 40) NAME 

0029 OPEN (UNIT-LUNRPT, NAME-NAME, TYPE-'NEW‘ ,REC0RDSI2E-l 33 , ERR-1000) 

0030 CALL BRMOD ( COE F , CR I T . WD , P LT , LAT ) 

0031 CLOSE (UNIT-LUNRPT) 

C 

0032 STOP 
C 

0033 1000 CONTINUE 

0034 WRITE (LUNTRM, 1020) 

0035 1020 FORMAT ( ‘0ERROR IN OPENING FILE’) 

0036 STOP 

0037 END 


C-2 




to 



*o 




to 

cs: 



Ui 

to 



to 

to 



O 

to 



C-* 

1 

to 


< 

in 

to 

■— 



Ui 




to 

Ml 

UI 


o 

to 

n 

CM 



V 


< 

ro 






K 




o 

-J 


UI 

Or: 

MJ 

_ 

Z 

z 


to 

< 


< 

8 

z 


—1 



M 

iA ^ 

V) Q 

lii ^ 

« to 

o to 


oc 

X 


o 

O 

< 


to 

to 


< 


wo 

tn to N 

(/) to «-• 

uj to 

to to -• 

- toto 


I I 
U3 ^ 


fsi rg 

« « 


CO lA ^ 

CO to 
UJ toto 
to toto 
Cl toto 

o I t 

< ro ro 


o 

pg 


M 


to 

.8 


u u o to 

-j _j _i tj 

z to to to 

o oo > 

o u u o 


to ^ o o o 

H 3r :> :r 

< to Qc to to 


lo to to i/> 
iri in 1/1 
n 


Z 

< 


to ^ 

to to CM 

to to 
ui to 
to to 

o toto 

“ I I 


to ^ 


o 

o 


< 

X 


^ -j 



to < « CM 


to 0 


lA tt 


to 

r- 1/1 ro — 

to 

CO to 

to 

1 •• 



5 ^ — to 

to 

to 

to 

CM Cf 



toto->to 

Ui 

to 

UI 

to ^ 



to to-*to 

to 

^ to 

to 

> 



totototo 

Q 

toto 

Q 


to 



o 

1 1 

o 

to 

z 

O 


< 

< 

^ 0 

< 

-t 

•— 

UJ 

UI K to 




to 

1- 

Z 

O < to 

UJ 


UJ 


•- h- to 


<* < 
to*- to 

»- < o 

toZ o 

O to to 

' . 0 a. 


to 

X 

3 


^ M 

* f! 
to 


« to 

< UJ H- 

« Z H Z 

to < O 


to 

z 

o ^ 

to n o 

z UJ 

UJ r> n ro 

X * *to * 

•— 4/1 m ro ro 

o ^ ^ ^ ^ 


to to to o 
oi ro ro cfi 

CM 

UI 

N 

•-* ^ to 

to 0 r«. in 

IM to to ^ 

ssss 

totototo 


to^toto 
to 0 ro VO 
to rsj toro 
toto -«to 

ssss 

• I I I 


to ' fSI 

< toto 
-I in — 


Of 

o 

o 

< 


to 

CO 


5 


to 

>• 

flC 

< 


^ ^ CM ^ 
« « « « 

oc Of «-af 


UI •- Z UJ to 

O0f<0 OD < 

UUZ> < ^ 


^ ro 
to ^ 
toto 
to to 
to to 

I I 


to 
. to 
to to 


toto 




I I 
m ro 


toCM 


3 

o 

Of 

to 

to 


< 

to 

z 

o 


o 

CJ 


to 

o 


§ 

to 


to 

to 

CM 


to 

■ 


< 

u 

o 


u 

< 


to 

< 


o 

♦- 


u 

to 

z 

o 


C-3 


Cl:*SV«:(32t,70)BRNAIN 


10:36:49 


02-OUN*81 


PAGE 1 


FORTRAN IV-PLUS V02-51E 
BRMOO . FTN /TR : BLOCKS/WR 


0001 


0002 

0003 

0004 

0005 

0006 
0007 


0008 

0009 

0010 
0011 
0012 
0013 


C 

SUBROUTINE BRMOD 

(COEF 

,CRIT,WD,PLT.LAT) 

c 

c 

SUBROUTINE BRMOD IS 

THE DRIVER FOR THE BASIC ROBERTSON 

c 

c 

CROP PHENOLOGY 

MODEL 

• 

c 

c 

r 

CALLING ARGUMENTS: 


c 

r 

NAME 

I/O 

MEANING 

L 

c 

COEF 

I 

COEFFICIENTS 

c 

CRIi 

I 

CRITICAL LEVELS 

c 

WD 

I 

WEATHER DATA 

c 

PLT 

I 

OULIAN PLANTING DATE 

c 

c 

LAT 

I 

LATITUDE 

c 

c 

c 

SUBROUTINES CALLED: 

DAYLG.ROB 

c 

c 

c 

CALLED BY: BRMAIN 


c 

c 

r 

LOCAL VARIABLES: 


L- 

c 

c 

NAME 


MEANING 

L 

c 

OS 

TRANSITION PERIOD 

c 

TMAX 

MAXIMUM TEMPERATUF. DEG. F 

c 

TMIN 

MINIMUM TEMPERATURE, DEG. F 

c 

PRECIP 

PRECIPITATION 

c 

DL 

DAY 

LENGTH 

c 

DIDO 

DAILY INCREMENT OF DEVEL.i OLD STAGE 

c 

OIDN 

DAILY INCREMENT OF DEVEL.; NEW STAGE 

c 

DIDI 

DAILY INCREMTNT OF DEVEL. { INTERPOLATED 

c 


FROM 

OLD AND NEW STAGE 

c 

c 

r 

CUM 

CUMULATIVE INCREMENT OF DEVELOPMENT 

c 

COMMON /LUN/ LUNCOE, 

LUNRES.LUNMET,LUNTRM.LUNRPT 


C 


INTEGER PIT 
REAL LAT 

DIMENSION COEF< 5,3,3 ),CRIT< 5,3 > ,W0( 3,366 > 

OS « 1 
CUM • 1.0 
C 
C 

DO 100 0D-PLT.366 
TMAX • W0< 1 ,0D) 

TMIN • WD<2,0D) 

PRECIP ■ WD( 3,00 > 

DL » DAYLG( LAT,0D) 

CALL ROB (TMIN, TMAX, OL,COEF,CR1T, OS, DIDO, DION, DIDl, CUM) 
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FORTRAN IV-PIUS V02-51E 
BRMOO.FTN /TR : BLOCKS/WR 


00U 

0015 

0016 
0017 


001 B 
0019 


WRITE (LUNRPT,50) 00 ,TMIN *T«AX ,OL ,0S ,DIDO,DIDN .DIOI ,CUM 
50 FORMAT ( * M 5 , 3F 1 2 . 4 , 1 6 , 4F 1 2 . 4 ) 

IF (CUM. CE. 6.0) RETURN 
100 CONTINUE 
C 
C 

RETURN 

END 
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